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ABSTRACT 


Additions  to  and  modifications  of  a  computer  program  for  the  numerical 
study  of  piezoelectric  surface  wave  propagation  are  described.  The  main 
addition  to  the  program  is  the  calculation  of  surface  wave  attenuation  using 
an  expression  based  on  a  perturbation  analysis  and  the  viscosity  tensor  which 
accounts  for  losses  due  to  the  Akheiser  mechanism.  Other  additions  and 
modifications  include  the  evaluation  of  stored  energy  expressions  and  changes 
in  the  format  and  normalization  of  the  output  data. 
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This  report  describes  additions  to  and  modifications  of  a  previously 
written  computer  program^  for  the  numerical  study  of  piezoelectric  surface 
wave  propagation.  The  extended  program  adds  the  :alculaticn  of  a  surface 
wave  attenuation  coefficient  o^,  representing  mechanical  losses,  for  a  general 
anisotropic,  piezoelectric  crystal  with  arbitrary  surface  cut  and  propagation 
direction  (including  the  loss  effects  in  the  metai  overlay  where  appropriate). 

Mechanical  loss  can  be  accounted  for  by  introducing  complex  elastic 
tensor  coefficients  c  at  the  outset.  The  resultant  complex  surface  wave 
velocity  would  thus  specify  the  attenuation  constant  as  well  as  the  phase 
constant.  This  procedure  is  computationally  disadvantageous  because  the  loss 
(imaginary)  terms  in  the  elastic  '’constants’'  are  absolutely  frequency -dependent, 
whereas  it  is  desirable  to  use  only  normalized  frequency  (xh)  when  there  is  a 
layer  of  thickness  h  contiguous  to  the  substrate  material. 

A  more  convenient  procedure  is  to  implement  a  perturbation  calculation 
in  which  the  lossless  (real)  propagation  velocity  is  first  found,  and  tl.-c 
attenuation  constant  a  is  then  computed  from  the  associated  strains  and  the 

P  (2) 

viscosity  tensor.  This  method  has  been  described  by  King  and  Sheard, v  '  and 

(3) 

includes  only  the  losses  due  to  the  Akheiser  mechanism'  (which  dominates 
the  temperature  dependent  losses  in  dielectrics). 

The  power  dissipated  per  unit  volume  by  the  wave  is  given  by  twice  the 
dissipation  function: 


2’*iJiUi  Vtsijsh 


where 


V2(\J+V  (2) 

are  the  tensor  strain  components  and  are  the  viscosity  tensor  components , 
The  viscosity  tensor  has  the  same  symmetry  as  the  elastic  tensor.  Note  that 
the  dissipation  function  depends  on  strain  time  derivatives  rather  than 
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particle  velocity  components  u.;  if  the  latter  entered,  a  translation  of  the 
crystal  would  result  in  loss . 

The  attenuation  constant  is  the  reciprocal  of  the  distance  (measured 
along  the  propagation  vector,  k)  over  which  the  intensity  of  the  wave  decreases 

-1  (4) 

by  a  factor  of  e  .  This  is  given  by'  1 


% 


2  ft  dx„ 
o 


(3) 


where  P  is  the  component  of  the  time  average  acoustic  flux  along  k  (the 
direction)  and  the  axis  is  normal  to  the  surface. 

The  computer  programs  of  reference  1  compute  P/a.  The  numerator, 
2j  i  dxQ,  is  computed  as  follows.  The  particle  displacement  components  of 
the  surface  wave,  in  one  medium  (substrate  or  layer)  take  the  form 


ui  =  ej«EB<q)6Wexp  -jf  T  ,®x, 
q  s  m=l 


m 


(4) 


In  (4),  y^  =  1 


and  Y2  =  0  are  direction  cosines  of  propagation  in  the 

^  is  related  to  the  q'th  transverse 


surface  plane,  and  are  independent  of  q .  y^ 
wave  number  by  y!  ^  =  -ja^.  The  use  of  this  y-notation  simplifies  the 

v* 

substitution  of  (4)  and  (‘  into  (1),  which  gives  the  result 


\  h  dx3  =  2 


1_ 

2v 


3 

E  E 

q,r  i,j,k, 


n,. 


ijkf. 


B(q,s;<rMq)vi(r)Q, 


qr 


a.(r)+0(q) 


(5) 


The  B's,  P's,  y's  and  a's  are  all  readily  available  in  the  computer  programs 
of  reference  1.  The  factor  Q  is  unity  in  a  crystalline  substrate  but  has  the 
value  =  exp  {uh/vs  To*^  +a^]}  -  1  in  a  layer  of  normalized  thickness 
ah.  q  and  r  are  summed  over  all  the  partial  waves  in  each  medium. 
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The  new  veision  of  the  program  evaluates  equation  (5)  as  well  as  P/u> 
and  then  calculates 


4  u&J  + 

4/+  dxs 

Luj  -J  J  layer 

Lu)  J  J  substrate 

[wl  substrate  1  ail 

layer 

The  terms  for  the  layer  are  simply  omitted  when  the  layer  is  absent.  Equation 

2 

(6)  displays  explicitly  the  x-  -dependence  of  the  attenuation  coefficient,  in  which 

2 

a  is  given  in  nepers/meter.  The  program  output  contains  1.715  a  /a-  which 
P  9  -  i  -  2  r.  _2  "  -4"  / 

gives  a  /1  in  dB  cm  Hz  .  It  also  prints  a  /i  (where  a  =  10  a  v  ,  v 
p  '  P  P  P  s  s 

in  m/sec),  thus  expressing  the  attenuation  in  dB  psec  Hz"2. 

Results  for  Y-cut  quartz  are  given  in  Figure  1.  Values  for  the  viscosity 
coefficients  were  taken  from  reference  5.  a'/f^  and  l/vg  =  S/2nf  are  plotted 
versus  direction  of  propagation  in  the  XZ  plane.  The  similar  shapes  of  the 
cui  ves  indicate  that  the  Q(=  3/ 2a)  has  relatively  small  anisotropy. 


30°  40°  50u  70°  SU°  90°  v- prop 
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APPENDIX  -  PROGRAMMING  EXTENSIONS 

The  attenuation  calculations  have  been  added  to  the  program^  for 
piezoelectric  crystals  in  the  presence  of  an  electric  conducting  plane  and  also 
to  the  cases  of  a  perfectly  conducting  isotropic  elastic  film  on  the  piezoelectric 
substrate.  The  corresponding  program  descriptions  of  reference  1  hold, 
subject  to  the  following  modifications. 

I .  Modifications  to  Input  Data  (see  $CQNST  data) 

A.  Changes  to  Existing  Input  Data 

The  elastic  constants  (C  )  must  be  input  in  the  following  order: 

W 

G  =  Cii»Ci2»C13’C14’C15’C16’C22>C23’C24’C25’C26’ 

C33*  C34*  C35*  C36’  C44’  C45‘  C46’  C55’  C56’  C66 

These  21  input  values,  together  with  the  21  transformed  elastic  constants 
(CA)  are  printed  out  in  the  order  shown  above  whenever  COEFF  =  .TRUE. 

B.  Additions  to  Input  Data 

In  order  to  compute  the  attenuation  coefficient  the  following 
parameters  must  be  added  to  the  input  data: 

SCONST-  Coefficients  of  Viscosity  (r.  ) 


(Optional)  Medium  A  (isotropic  layer)  -  input  two  values  of  rj 
ETAA  = 


which  are  used  to  generate  a  twenty-one  (21)  element  array  whose 
elements,  in  standard  reduced- subscript  notation  are: 
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V 


II "  11 

12  =  12 
"13  s  12 


'14  =  ° 


=  0 
Id 


16-° 

22  =  11 


23  12 


"24  =  ° 


25  =  0 


26 


=  0 


r33  =  "ll 


34 


=  0 


r,-  =  0 
3d 


'36 


=  0 


44  "  44 

45  =  ° 

'46  =  0 
55  ~  44 


56  =  0 


66  ~  44 


—  "  -  9  r 

12  11  44 


i 


f 


* 


« 


where 
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Medium  B  (substrate)  -  input  all  21  values  of  n  in  the  order 
listed  above 


ETAB  -  rjn,  m»2>  t113,  •  •  ^22,r'23’  *  *  **  ni 


66 


Wherever  printout  is  requested,  both  the  input  constants  and 
the  transformed  constants  are  tabulated  in  the  order  shown  above. 


$INPUT 


Input  Name  Equation  Name  Definition 

FATEN  --  A  logical  parameter  which  controls  the 

calculation  of  the  Attenuation  Coefficient 

•TRUE.  -  Calculate  the  Attenuation 
Coefficient  (input  '"‘s  in 
the  $CONST  data) 

.FALSE .  -  Do  not  calculate  the 

Attenuation  Coefficient. 

(Nominal  Value  -  .FALSE  ) 


II .  Modifications  to  Output  Data 

A.  Changes 

1.  The  input  elastic,  piezoelectric  and  dielectric  constants 
are  included  in  the  printout  of  the  transformed  constants . 
(The  elastic  constants  have  been  re-ordered  for  both 
input  and  output  —  see  I. A.) 

2.  The  FINAL  ANSWERS  are  now  normalized  as  follows: 

All  field  quantities  (  u,  T/oj,  S/w,  E/ai,  D/o)  are  first 
found,  as  before,  with  arbitrary  normalization .  From 
these  amplitudes  a  value  of  (P^/tl)  is  found.  The  field 
quantities  which  are  printed  out  are  normalized  through 
division  by  the  factor 
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t- 
■  & 
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> 


Medium  8  (substrata  —  input  ail  21  values  of  tj  in  the  order 
listed  above 


ETAS  nir  T*i2’  ^13'  *  V,22*  *  "  *a  ^ 
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Wherever  printout  is  requested,  both  the  input  constants  anti 
the  transformed  constants  are  tabulated  in  the  order  shown  above. 


SjgjPUT 


Input  Name  Equation  Name  Definition 

FATEN  --  A  logical  parameter  which  controls  the 

calculation  of  the  Attenuation  Coefficient 

-TRUE-  -  Calculate  the  Attenuation 
Coefficient  (input  r)*&  in 
the  $CONST  data) 

-FALSE.  -  Do  not  calculate  the 

Attenuation  Coefficient. 

(Nomiual  Value  =  -FALSE.) 


II.  Modifications  to  Output  Data 
A.  Changes 

1.  The  input  elastic,  piezoelectric  and  dielectric  constants 
are  included  in  the  printout  of  the  transformed  constants. 
(The  elastic  constants  have  been  re-ordered  for  both 
input  and  output  —  see  I. A.) 

2-  The  FINAL  ANSWERS  are  now  normalized  as  lullows: 

All  field  quantities  (  u,  T/iju,  S/uj,  E/uj,  D/uo)  are  first 
found,  as  before,  with  arbitrary  normalization.  From 
these  amplitudes  a  value  of  (P^/u)  is  founc.  The  field 
quantities  which  are  printed  out  are  normalized  through 
division  by  the  factor 


s 


Tree  Sms r  lime  average  power  Oow  components,  (?j/») 
3c-d  (P  ^ :.}  for  media  A  and  B,  are  each  normalized  lo 
Re(F./j;)g.  Wills  the  aiwve  corscalization  (P . / lakes 
ik  value  unity . 


B.  Adatioas  (included  in  FINAL  ANSWERS  only) 

1-  In  the  substrate  te  normalized  transverse  wavenumbers 

(/ >  /A 

r  '  and  the  corresponding  amplitudes  ?;  ’  for  each 
partial  wave  are  printed. 

2.  A.  tabular  summary  of  the  following  quantities  is  printed 
for  each  value  of  propagation  -iirection  angle  (-)  when 
A  v?£  0 


* 


i 


v  I/v 

s  s 

(degrees)  fns/sec)  (sec /in) 


^  /  * 

P  ?  ?  ? 

(dB/cm  Hzz)  (dB/psec  ilz“) 


-1 

tan  x 


tan  *x  is  the  angle  between  the  i:  vector  and  the 
mechanical  power  flow;  thus 


Re 

l“fj 

A+Re 

[!i 

[TjB  layer 

Re 

Pll 

i  J  ] 

A+Re 

TT”  present 

T/b 

Re  j 

|P21 

ijJ 

B 

i 

Re  j 

lPll 

1  '! 

i 

|b 

i 

V. 


layer 

absent 
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For  lbs  piezoelectric  substrate  with  a  perfectly 
conducting  plane  at  a  distance  rn  =  0  and  at  a  distance 
jJi  =  *  the  toll  owing  stored  energy  terms  are  printed: 
?£,  P  P-y  P^,  These  terms  are  defined  and 
explained  in  the  following. 

When  the  constitutive  relations 


+-+■  ■♦i 

S  =  s  :  T  +d 


~~t~  ♦♦  +♦•  — ♦ 

D  =  d  *  T  +  e  •  E 


are  used,  the  total  stored  energy  is  given  by  the  sum, 
?E+  2PM+PD*  where 


*  H-iffr  •***:?.  t) 

PD  =ifr 


dv 


•  s  •  E  dv 


The  appropriate  computational  expressions  are: 
a)  Elastic  energy  term  evaluated  using  stresses 
6  4 


Bw  B*W 

-if  f  fWf^s 
E  2  u=i  k,%i  1  J  *1  a<  +a° 


b)  Mutual  energy  term 


Jn  ti  ci  U\  R  ^ 

2P  =  Re  >  f  \  V  fW  c  W  ^  £ 2 - 

m  A*  A  p  p.  i  sm  mi  (-C)  ,  *(k) 

i=l  m=l  k=l  £=  1  a  +  a  '  ' 
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c)  Electric  energy  term,  dielectric  constant  at 
constant  stress 


1  J 

*D-i  £  , 

m,  n=l 


U)  *  (k)  „T 

^n  'mn  (l)  ,  “(k) 
o  +  a  v  ’ 


+ 


ih  =  0 

jJi  =  ^ 


,  Xt)  (k) 

where  f>  \  g^'. 


and  h: 


(<> 


are  given  below . 


An  alternative  set  of  constitutive  relations  is 


-*->  -H^r  — y 

T  -  C  :  S  -  e  •  E 


— v  -+-tr  +r  — T 

D  =  e:  S  +  •  6 


When  these  are  used,  the  total  stored  energy  in  the 
piezoelectric  substrate  decomposes  into  +  ?£, 
where 

1  fr**  -f-y 

Pi=2  J1  *C  *S  dv 

and 


dv 


The  appropriate  computations  expressions  are: 
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d)  Elastic  energy  term  evaluated  using  strains 


6  4 


■,-4  E  f 

1  2  i,j=i  1  j  « 


(i)  *{k)  B(V<k> 

1  J  «  aW+a*(k> 


e)  Electric  energy  term,  dielectric  constant 


at  constant  strain 


p, = |  E  E  «(i)  g*(lt)  6s  -(t>  B*w 

2  2  m,n=l  kfUl  m  “  mtl  aTO  +  0*W 


ah  =  0 


.  „  B<t>  ^  I2  ’  *- 


For  the  above  five  energy  quantities  we  define: 
C.j  =  elastic  constants  (program  input) 

s„  =  (C'Vj 


e. -s.. 
ij  jk 


(e„  =  program  input) 


=  ei;  + 


^  dik  ejk  ^  = 


program  input) 


while  the  f^,  gj^,  and  h^  are  defined  by: 


i  '  °m 


^■l,it,c-jc11-«wc15]+4<)c-jc16-«wc14] 


+  S3>C-jCl5-o<%3]-,  ^C-i  cu  -  aW  e3i] 
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